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ABSTRACT 

We examine the possibility of constraining the cosmological mean mass and dark en- 
ergy densities by an application of the Alcock-Paczyhski test on redshifted 21-cm maps 
of the epoch of reionization. The 21-cm data will be provided as a function of frequency 
and angular positions on the sky. The ratio of the frequency to angular distance scales 
can be determined by inspecting the anisotropy pattern of the correlation function 
of brightness temperature. We assess the sensitivity of the distance ratio to the cos- 
mological parameters and present a technique for disentangling geometric distortions 
from redshift space distortions caused by peculiar motions. 

Key words: cosmology: theory - intergalactic medium -large-scale structure of the 
universe 



1 INTRODUCTION 



Observations of the cosmic microwave background (CMB) 
and Type la Supernova place significant constraints on the 
physical parameters governing the evolution of the cosmo- 
logical background (Spergel et al. 2003; Knop et al. 2003; 
Riess et al. 2004). The main parameters constrained by 
these observational data are the mean matter (mass) den- 
sity, f2 m , and the density of negative pressure energy (dark 
energy), Q v . Perhaps the most intriguing implication of the 
observations is the need for a dark energy component which 
dominates the energy budget in the Universe today (Ratra 
& Peebles 1988; Wetterich 1995; Coble, Dodelson & Frie- 
man 1997; Turner & White 1997; Wang et al. 2000; Cald- 
well & Doran 2004; Kunz et al. 2004; Caresia, Matarrese 
& Moscardini 2004). By and large the constraints derived 
from these observations are sustained by other estimates 
based on analyses of cosmic flows (Strauss & Willick 1995; 
Nusser, Willick & Davis 1997; Zaroubi et al. 2001), abun- 
dance of rich galaxy clusters and its evolution (Bahcall et 
al. 2003; Ikebe et al. 2002; but see Vauclair et al. 2003) the 
mass-density power spectrum derived from galaxy redshift 
surveys (Percival et al. 2002; Zehavi et al. 2002; Tegmark 
et al. 2004), and the Ly-a forest (Croft et al. 2002, McDon- 
ald et al. 2000, Nusser & Haehnelt 2000, McDonald et al. 
2004, Viel, Weller & Haehnelt 2004). Nevertheless all these 
constraints are derived from measurements at either very 
high (the CMB) or low and intermediate redshifts (SN and 
other data), excluding the wide redshift range the last scat- 
tering surface of CMB photons until z ~ 4. Here we present 



in detail a technique for constraining the cosmological pa- 
rameters by an application of the Alcock-Paczyhski (AP) 
test (Alcock & Paczyhsky 1979; see also Hui, Stebbins & 
Buries 1999; Ballinger, Peacock & Heavens 96, da Angela, 
Outram & Shanks 2005) on future data at z ~ 6 — 30. These 
data are maps of the redshifted 21-cm emission/absorption 
lines produced by H I in the high redshift universe (Field 
1959; Sunyaev & Zel'dovich 1975; Hogan & Rees 1979; Sub- 
ramanian & Padmanabhan 1993; Madau, Meiksin & Rees 
1997). The 21-cm line is produced in the transition between 
the triplet and singlet sublevels of the hyperfinestructure of 
the ground level of neutral hydrogen atoms. A patch of H I 
would be visible against CMB when its spin temperature 
T s differs the CMB temperature, T CBR . Various mechanisms 
exist for raising T s significantly above T CBR during the era 
of reionization. As a result a significant cosmological signal 
from the era of reionization should be expected (Field 1958 
& 1959; Scott & Rees 1990; Madau, Meiksin & Rees 1997; 
Baltz, Gnedin & Silk 1998; Tozzi et al. 2000; Ciardi & Madau 
2003; Chen & Miralda-Escude 2004; Gnedin & Shaver 2004; 
Nusser 2005; Ricotti, Ostriker & Gnedin 2004). 

We advocate the the application of the AP test on the 
correlations measured in three dimensional maps of 21-cm 
emission. These correlations are readily expressed in terms of 
angular separations and frequency intervals. In order to de- 
rive the correlations in terms of real space separations (mea- 
sured for example in h _1 Mpc) one needs to know the mean 
density parameter, fi m , and the dark energy density Q v and 
its temporal evolution (or equivalently its equation of state). 
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Constraints on these parameters can then be obtained by 
demanding isotropy of correlations in real space. There are 
two main obstacles in implementing this program on future 
data. The first is foreground contamination (e.g. Oh & Mack 
2003; Di Matteo, Ciardi & Miniati 2004) which we will not 
discuss here. The second is redshift distortions. The redshift 
coordinate of a patch of matter is given by the sum of its 
true real space coordinate and the its line of sight peculiar 
velocity (deviation from Hubble flow). The distribution of 
matter in redshift space should therefore be different from 
the real space distribution. According to the standard cos- 
mological paradigm, the observed large scale structure have 
formed by gravitational amplification of tiny initial fluctua- 
tions (Peebes 1980, Peacock 1999). In late time linear the- 
ory the gravitational growth of mass density fluctuations is 
inevitable associated with coherent peculiar motions. There- 
fore, over large scales where linear theory is generally valid 
the density in redshift space appears more clustered than 
in real space. Further, because because redshift and real 
space coordinates differ only in the line of sight direction, 
the structure in redshift space would appear anisotropic. On 
small scales random motions are important, causing a sup- 
pression of clustering in redshift space. Here we restrict the 
analysis to linear theory which is valid at the high redshifts 
(z ~ 10 — 30) and the physical scales ( > 1 h _1 Mpc) we con- 
sider. Redshift distortions, in the linear as well the non-linear 
regime, have been investigated in great detail in relation to 
the galaxy redshift surves (Peebles 1980; Davis & Peebles 
1983; Kaiser 1987; Hamilton 1992; Cole, Fisher & Weinberg 
1994; Nusser & Davis 1994, Fisher & Nusser 1996; Taylor 
& Hamilton 1996; Zaroubi & Hoffman 1996; Hatton & Cole 
1998; Desjacques et al. 2004) and to a lesser extent in rela- 
tion to 21-cm maps (Illiev et al. 2002; Furlanetto, Sokasian & 
Hernquist 2004; Barkana & Loeb 2004, Bharadwaj & Saiyed 
2004)). 

In this paper we show that the pattern of redshift space 
anisotropics is quite distinct from the geometric distortions 
introduced by an erroneous choice of the cosmological pa- 
rameters. Further, we show that an application of the AP 
test will provide useful constraints on the cosmological pa- 
rameters if the correlation function of the cosmological 21- 
cm brightness temperature could be estimated to within a 
20% accuracy. 

The outline of the paper is as follows. In §2 the relation 
between 21-cm emission/absorption and the gas density is 
briefly summarized. §3 describes in some detail redshift dis- 
tortions in the linear regime and the expected anisotropy 
the density correlation function in redshift space. In §4 the 
dependence of geometric distortions on the cosmological pa- 
rameters is explored. §5 shows how to distinguish between 
geometric and redshift space distortions. In §6 we discuss 
additional application of the AP test with 21 cm maps. We 
conclude with a short discussion in §7. 



(Wild 1952). The differential brightness temperature (here- 
after, DBT) against the CMB of a small patch of gas with 
spin temperature T s at redshift z is (e.g. Ciardi & Madau 
2003), 



8T h = 16mKi 



(1 + 5) (l- 
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where the fraction of H I in the patch. The quan- 

tity 8 — p/p — 1 is the density contrast of the gas. The 
fraction x HI may depend on position, r, and the density con- 
trast, 5 (e.g. Knox et al. 1998, Buscoli et al. 2000, Benson et 
al. 2001). This dependence is mainly dictated by the nature 
of the ionizing radiation and the spatial distribution of its 
sources. If the early staged of reionization are dominated by 
X-ray photons which have very large mean free path then 
x HI at any point in space is determined by the equilibrium 
between photoionization and recombinations (e.g. Ricotti & 
Ostriker 2004). In contrast, UV photons from stellar sources 
have a short mean free path. Therefore, in the presence of 
UV photons alone rionization proceeds in a patchy way with 
x m very close to zero in ionized regions and unity otherwise. 
We model reionization in the presence of X-rays and UV 
photons as follows. We assume that the explicit dependence 
of x HI on 8 outside the regions ionized by UV is linear so 
that, 



x H1 (r,8) oc M(r) 



1 + 



d In x E 
~~£)8~ 



(2) 



where M(r) is termed the mask and is either zero or unity. 
The expansion of x HI to first order in 8 is justified here since 
we only consider correlations on scales where the density 
contrast is small. 

In §5 we will present a formalism for an application 
of the AP test using the correlations of 8T b . In develop- 
ing that formalism we will assume that the mask structure 
(correlation) function is insensitive to peculiar velocities and 
also uncorrelated with the 21 cm signal. We will see that 
these conditions are valid in the limit of large filling fac- 
tors (F{ > 0.5) of ionized regions At the early stages of re- 
ionization, when the filling factor is small enough, it will be 
possible to identify large unionized regions in future maps 
of 21 cm (see figure 3 of Ciardi & Madau 2003). At these 
early stages, all complications related to the mask, M, can 
completely be avoided by restricting the analysis to these 
regions. In an x-ray pre-reionization scenario, X-rays from 
accreting black holes contribute to partial ionization before 
the appearance of significant UV sources (Ricotti & Ostriker 
2004). During this phase, the mask is unity except in the im- 
mediate vicinities of the X-ray emitting sources (Zaroubi & 
Silk 2005) . The number density of X-ray sources is expected 
to be low (Ricotti & Ostriker 2004) so that role of the mask 
is negligible in this case. 



2 THE BRIGHTNESS TEMPERATURE 

Intensities, I{v), at radio frequency are expressed in terms 
of brightness temperature defined as T b = I(y)c 2 /2ku 2 , 
where c is the speed of light and k is Boltzman's constant 



3 REDSHIFT DISTORTIONS 

For simplicity of discussion we consider a region in the shape 
of a cubic box at redshift z. Further, we assume the "distant 
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observer limit" (e.g. Kaiser 1987, Zaroubi & Hoffman 1996) 
according to which the box is small compared to its distance 
from the observer. We work with a Cartesian coordinate 
system defined by the axes of the box assuming that one of 
these axes is in the direction of the line of sight to the box. 
The box is assumed to be comoving with the Hubble flow 
and large enough so that the center of mass of matter inside 
it is also comoving with the Hubble flow. 

In the coordinate system attached to the box, let v and 
r be, repsectively, the physical peculiar velocity (deviations 
from Hubble flow) and physical real space coordinate (i.e. 
Eulerian) of a particle, both expressed in km s~ . The red- 
shift space coordinate is defined as s = r + vyl, where 1 is a 
unit vector in the line of sight direction and v« = (1 • v) is 
the line of sight peculiar velocity*. It convenient to denote 
the components of r parallel and perpendicular to the line 
of sight by ry and rx, respectively. Using similar notation 
for s we find si = rx and sy = ry + wy. 

Only in the absence of peculiar motions (both thermal 
and cosmological) the density contrast, S, appearing in (1) 
refers to the actual real space density constrast, S T . In the 
presence of peculiar motions 5 is equal to the gas density in 
redshift space, S s . This is given by 

l + «f(sx,sy) = 

J dry [l + 5 r (rx,rn)] G [ry + U(r ± , ry) - *„] , (3) 

where we have decomposed the line of sight peculiar veloc- 
ity into a smooth component U(r) and a random thermal 
component, f t h having a gaussian probability distribution 
function, G(vth), of mean zero and rms oth- In the following 
we assume that cr t h is negligible so that G reduces to a Dirac 
delta function. In this case we obtain 



I -,,» = [I f ,,'(r)j ( l-g!L 



where ry is determined by ry + U(r±,r\\) = sy 
(4) to first order in 5 and U yields 



(4) 



Expanding 



8 s W 



8 ( r ) ~~ -T~- 

dm 



(5) 



Note that to first order this last relation could equally be 
expressed as a function of s rather than r (Nusser & Davis 
1994). The term dU/drn introduces anisotropies in the den- 
sity field in redshift space. Supplemented with a relation 
between the velocity and density fields of the gas, the rela- 
tion (5) allows us to determine S s from S T . On large scale, a 
density- velocity relation between the real space mass density 
and velocity fields can be obtained using linear perturbation 
theory (e.g. Peebles 1980; Peacock 1999) . This relation reads 



(6) 



As an illustration of redshift distortions we show in figure 
(2) slices of the density fields in real and redshift spaces in 

* The observed redshift coordinate also contains contributions 
from the cosmological expansion and from the prcculiar motion 
of the observer. However, these contributions do not amount to 
any systematic effect. 



the bottom left and top left columns, respectively. The real 
space density field is a gaussian random field generated with 
the linear density power spectrum of the cold dark matter 
cosmology (CDM) with f2 m = 0.3. The linear relations (5) 
and (6) have been used to derive S s for a line of sight is 
in the y-direction. The density fields are shown in a slice of 
128 h _1 Mpc on the side. To improve the visual presentations 
the fields have been smoothed with a gaussian window of 
3 h _1 Mpc width. The normalization is arbitrary and the 
contour spacing is 2. A comparison between the top left and 
bottom left panels reveals significant differences between S B 
and 5 T . The lack of isotropy in the contour maps of 5 s is also 
clear (see Kaiser 1987 for more details). 

The relation (6) can easily be generalized to account 
for any bias between S T and the actual mass density, 5^- 
To first order in the density contrast such a bias can be 
described by S T = b5 T m where & is a constant. Then the re- 
lation (6) is modified by multiplying the right hand side by 
b. Hamilton (1992) used the linear relations (5) and (6) in 
order to express the redshift space density correlation func- 
tion £ s =< <5 s (ro)5 s (r + r > r o in terms of the real space 
correlation £ r =< S T (ro)S T (ro + r > r o as follows, 



e s (m 



£o(r)P,(/x) + Hr)P2(») + &(r)P 4 (/*) 



(7) 



where ry and fx are the components of r parallel and per- 
pendicular to the line of sight, /1 = ry/r is cosine the an- 
gle between r and the line of sight, and p is the I th or- 
der Legender polynomial: Po = 1, P2 = (SfJ 2 — l)/2, P4 = 
(35^ 4 -30^ 2 +3)/8, andP 6 = (231^ 6 -316^ 4 + 105^ 2 -5)/16 
(we list P5 for subsequent use). Further, 



6M = (3/3+ 7 /3 3 j [rw-e] 



(8) 
(9) 
(10) 



where (3 = ^ /b, 



3r" 



and f r 



4 GEOMETRIC DISTORTIONS 

Twenty one centimeter maps will be readily provided in 
terms of the frequency along the line of sight direction and 
angular positions on the sky in the perpendicular plane. An 
observed frequency interval Au corresponds to a physical 
separation (in km s _1 ), Ary, of 



Arii = 



1420Hz 1 + 2 



(11) 



This relation is independent of the cosmological parame- 
ters. The projected separation (also in km s _1 ), Arx, cor- 
responding to an angular separation of A# is 



Arx = AdHD A {z) 



(12) 



where Da = DA(fim, O v , z) is the angular diameter dis- 
tance and H(Q m ,Q v , z) is the Hubble function evaluated 
at redshift 2. Therefore, for a given angular separation one 
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needs to know the cosmological parameters in order to de- 
termine Arx- The dependence HD\ on the cosmological 
parameters is illustrated in figure (1) which shows the quan- 
tity 1/(HDa) for cosmological models with and without a 
dark energy component (Lima & Alcaniz 2000). We assume 
a dark energy equation of state of the form P v = wp v , where 
P v and p v are, respectively, the pressure and density of the 
dark energy component, and w is constant. This equation 
of state implies p v oc (1 + 2 ) 3 ( 1+1 ") and so w = — 1 describes 
a cosmological constant (constant p v ). We do not consider 
here w < — 1. The figure shows 1/(HDa) normalized to its 
value^ for fi m = 0.3, Q. v = 0.7 and w = — 1. The solid line 
in the figure is for an open universe universe with £l v = 0. 
All remaining curves are for a fiat universe (O m + Q v — 1) 
and they correspond to several values of w as indicated in 
the figure. The distance ratio is clearly sensitive to the as- 
sumed cosmological parameters. For the models plotted in 
the figure the maximal difference is about 40% at fi m = 0.2. 
The panels to the right of figure (2) illustrate the nature of 
the geometric distortions. The top (bottom) panel show the 
density contours obtained by stretching the vertical (hori- 
zontal) axis of the slice in the left bottom panel 25%. The 
distinction between geometric and redshift space distortions 
is readily seen by comparing these maps the density field 
in redshift space shown in the left top panel. Redshift dis- 
tortions amplify the density fluctuations, while geometric 
distortions merely stretch the contours. 

In the absence of peculiar motions and foreground con- 
tamination the correlation function of the brightness tem- 
perature is isotropic for the true f2 m , fi v , and w. However, we 
expect redshift distortions to be present as density fluctua- 
tions are associated with peculiar motions under the action 
of gravity. We show in the next section how to disentan- 
gle the geometric from redshift distortions, allowing an es- 
timation of the cosmological parameters from the isotropy 
pattern of the correlation function. 

In the remainder of this section we derive several rela- 
tions for later use. The relation between redshift and real 
space correlations is invariant under an isotropic scaling 
transformation of the form r — > gr, where g is a constant. 
Therefore, without loss of generality the effect of an incor- 
rect choice of the cosmological parameters can be described 
by a transformation of the type f || = r||/(l+a) and fx = rx, 
where (1 + a) is the ratio of the true value of HD\ to the 
assumed value. For clarity of notation, we place a tilde over 
functions defined in the space of fy and fx - If / = / (ry , rx) 
then 



/(fy,fx) = /[(ry,rx)] , 



(13) 



We can easily verify that the correlation function of d com- 
puted in the space of (f||,fx) is 

f(fy,rx)=f[(l+a)fy,fx] • (14) 
f(r) is independent of the direction 



According to (13), if / 
of r, then 



/( 



f(r) 



(15) 



t In this paper the numerical values of Q m and Sl v always corre- 
spond to z = 0. 




Figure 1. The dependence of the distance ratio 1/(HD\) at 
z = 20 on the assumed cosmological parameters. The ratio is 
normalized to its value at Q m = 0.3, Q v = 0.7, and w = — 1. 
The solid curve is for an open universe with matter content only. 
All other curves correspond to universes with f2 m + Q v = 1 for 
various values of w, as indicated in the figure. 



where r — ^J(l + a) 2 r jj + f \ = r^J (1 + a) 2 fi 2 + (1 
To first order in a we find 



/(n|,fx) = f(r) 



, 2 din /(f) 
1 + ap? J w 



dlnf 



A 2 ) 



(16) 



(17) 



We will also need the first order (in a) relation 
fj 2 = jl 2 + 2a.fi 2 (1 — ft 2 ) . 
between jl = ru/f and p — r\\jr. 



5 DISTORTIONS IN TWENTY ONE CM MAPS 



We assume that T s T CBR so that STh oc x H1 (1+<5 S ) up to a 
multiplicative constant. We further substitute the expansion 
(2) for x m and write a; HI (l + <5 s ) cx M(r) [1 + (6 - l)S r + S s ] 
to first order in S, where 6 = 1 + d\nx HI /d6. Note that 
x HI is determined by the real space density <5 r , hence the 
appearance of 8 1 . For simplicity we will assume that the 
data yields the quantity 



d(r) =Af(r)(l + A s ) , 



(18) 



where A s = (b—l)S I +8 s . Using the linear relation (5) we find 
that A s = (b-l)<5 r +<f = bS'-dU/dr which can be identified 
as the redshift space density contrast corresponding to the 
real space density field A r = b8 r . 

We write the data-data redshift space correlation func- 
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Figure 2. An illustration of the redshift and geometric distor- 
tions. The bottom left panel with the isotropic contours shows 
the density field in real space. The top left is the corresponding 
redshift space density field obtained from the linear relations (5) 
and (6) for a line of sight along the vertical axis. The top right and 
the bottom right panels illustrate geometric distortions parallel 
and perpendicular the line of sight, respectively. These geometric 
distortions are for a 25% difference between the assumed and true 
distances 25%. 




0.5 1 1.5 

Log 10 [r h/Mpc] 

Figure 3. The correlation tp as a function of separation, for two 
values of the filling factor of ionized regions, as indicated in the 
figure. The upper (Ff = 0.13) and lower sets (F{ = 0.36) of curves 
are for sources identified as 4<r and and 3.5<r peaks in the density 
field, respectively. The solid lines show the correlation when the 
mask of ionized regions is placed around the peaks, while the 
dotted lines are for masks displaced away from the peaks. 



tion as 



^ B =<dd>=C M (r)[l+e(r)] 



(19) 



where C M =< MM > is the structure (correlation) func- 
tion of the mask, and £ s is now the correlation function of 
A s . We have also assumed that the mask x H1 and 5 are un- 
corrected. This assumption is justified once the filling factor 
of ionized regions is larger than nR^ where n is the number 
density of sources and R c ~ h _1 Mpc (comoving) is the co- 
herence length of the correlations between the sources and 
the gas. This condition is easily satisfied as a single source is 
capable of ionizing relatively large regions around it. To test 
this assumption we have resorted to our random gaussian 
density field having a ACDM power spectrum (see figure 2). 
We identify the sources as peaks in the density smoothed 
with 1 h _1 Mpc and assume that each source ionizes a region 
of a given radius around it. These regions define our mask. 
A comparison between the correlation of the mask times the 
density field when the mask is placed around the sources and 
when it is placed randomly in the box serves as an indica- 
tion to the the degree of correlation between the mask and 
the density field. In figure (3) the solid lines are the corre- 
lation function %j) (in real space) for a mask placed around 
the sources while the dotted are the correlations when the 
mask is placed randomly in the box. The upper and lower 
sets of curves correspond to two values of the filling factor, 
Ff, of ionized regions. The values Ff = 0.13 and 0.36, re- 
spectively, correspond to sources identified as 4a and 3.5<r 



peaks. The differences between the dotted and solid curves 
for each value of the Ff are negligible. 

We will write now the expressions describing the 
anisotropy of the data-data correlation function in the pres- 
ence of geometric and redshift distortions. We restrict the 
analysis to the linear regime 5 <C 1 and to first order in a. 

We write | s (f||,f_L) = £ s ((l + a)f||,r±) and use (7) 
to express £ s in terms of the moments £;(r) with (r = 

(1 + a) 2 f jj + and Legendre polynomials evaluated at 

fl = (1 + a)rn/r. Then we use (16) and (17) to expand 
and Pi as functions of f and ft to first order in a. The final 
result of this straightforward algebra is 



i> s (r h r±) = ^ ^i{r)Pi{ft) 



!=0,2,4 



dlnV>;(r)_ 2 . ~2,dP l (ft) 



-ft, Piifi) + ft{l -ft )- 



din? n ^ ^ ^ ' dp. 
= [Mr) + aai(r)]Pi(ft,)+aa 6 (r)P 6 {ft) (20) 
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where ip = C M (1 + Co) an d ipi = C M £, t f° r I > 0, and 

1 , dlnV> n , (2 2 dlnV>,\ 

2 , dln^n , /2 11 dlnV 2 \ 
, /20 4 din </>. \ 

, /24 12dlnV 2 \ , /20 39dhn/> 4 \ 
/40 lOdJn^N 

The visual impression between geometric and redshift 
distortions seen in figure (2) is quantified by the relation 
(20). The moments ai and ipi, respectively, characterize the 
geometric and redshift distortions. A term containing Pq is 
missing from the redshift distortions of the correlations. One 
could determine the correct value of the distance ratio, a, 
by minimizing the square of the moment corresponding to 
P6, independent of the details of the mask structure function 
C M . As we shall see below other moments can also be used to 
infer a. Let us inspect the moments a; and tpi. The moments 
depend on the linear matter correlation function which we 
obtain from the power spectrum of the ACDM model with 
flm = 0.3. Figure (4) shows the ratio of the moments cii/ipi 
for / = 0, 2, and 4. Also shown is multiplied by a factor of 
10 for the sake of clarity. To illustrate the differences arising 
from geometric and redshift distrotions, we assume that the 
mask structure function, C M , is unity. The ratio |a /V> is 
small except at separations of r ~ 60 h _1 Mpc where mea- 
suring the correlations is probably hard. The moments ai 
and <24 are larger than ip2 and ip4 by a factor of 2 and 10, 
respectively, making the anisotropy pattern a very sensitive 
function of a. 

All moments J i/> s (/i)F;(/i)d/i of the temperature corre- 
lation are expressed in terms of the functions Cm, £ , the 
redshift distortion parameter /3 needed to express £, in terms 
of £ , and the geometric distortion parameter, a. A compar- 
ison between the moments should allow us to constrain a, 
the bias parameter /3, and the form of the mask. As ex- 
plained previously the moment corresponding to arises 
only because of geometric distortions. Therefore, a can be 
constrained by minimizing the square of this moment inde- 
pendent of the mask and redshift space distortions. Further, 
one can assume parametric forms for Cm and £ . The param- 
eters of these forms as well as a and f3 can then be obtained 
by fitting the moments as a function of r. 



5.1 Validity of the approach 

During the early stages of reionization the filling factor, 
Ff, of ionized regions is small and the effect of the mask 
on the temperature correlations is negligible. Further, one 
can restrict the correlation analysis to large neutral patches 
extracted from 21 cm maps (figure 3 in Ciardi & Madau 
2003). In the intermediate and late stages of reionization 
(Ff > 0.5), it will be increasingly difficult to extract large 
neutral regions from the data (Nusser et. al. 2002). During 



these stages, the temperature correlation from the full maps 
will probably need to be considered. The method outlined 
above is capable of treating these stages where the mask is 
to be considered. In developing the method we have made 
two important assumptions related to the mask. The first 
is that the mask-temperature correlation is negligible. The 
second is that peculiar motions do not introduce any sys- 
tematic anisotropy in the shape of the mask boundaries, 
i.e., the ionization fronts. Once the typical size of ionized re- 
gions becomes larger than the correlation (coherence) length 
( < 1 h~ Mpc at z ~ 10) of mass distribution the mask- 
temperature correlation is expected to be negligible. Figure 
3 illustrates that the mask-signal correlation is indeed neg- 
ligible once the ionized volume is large enough (correspond- 
ing to Ff > 0.13 in the example given in the figure). Let us 
now justify neglecting the effect of peculiar motions on the 
mask boundaries. The relative displacement of the bound- 
ary is v/(aHR) where R is the typical comoving size of the 
ionization front, and v is the typical difference in physical 
peculiar velocity between two points at which a line of sight 
intersects the boundary. Also, a = 1/(1 + z) and H is the 
Hubble function at redshift z. For simplicity consider Q — 1. 
In this case v = voa 2 H / Ho where vo and Ho are v and H at 
z = 0. The relative displacement is then avo/(HoR). Taking 
HoR = lOOOkm/s which is appropriate for Ff > 0.5, a — 0.1, 
and vo = lOOkm/s one gets a relative displacement of a few 
per cents. There is also an important factor which reduces 
any systematic distortions of the ionization fronts. For large 
ionized regions the the relative velocity v can equally be 
positive as much as negative. This is because UV ionized 
bubbles are significantly larger than the coherence length 
(~ 1 h _1 Mpc at z ~ 10) of the dense regions near which the 
sources are likely to reside. Because of this symmetry in v, 
the mean distortions in the boundaries of bubbles tend to 
cancel out. 



6 ALTERNATIVE APPLICATIONS OF THE 
AP TEST 

So far we have considered the anisotropies of the correla- 
tions in the diffuse neutral component. We now consider 
two alternative applications which are also promising and 
are practically unaffected by foreground contamination. 



6.1 Dense peaks and minihalos 

Consider a UV dominated reionization scenario. At any red- 
shift, most of the 21 cm signal comes from the regions that 
have not yet been ionized. However, the already ionized re- 
gions are likely to contain high density peaks in which the 
neutral hydrogen is dense enough to produce significant 21 
cm emission. For sufficient instrumental spatial resolution 
these peaks would show localized 21 cm emission engulfed 
in regions with practically no 21 cm signature at all. These 
peaks are particularly useful in the absence of 21 cm sig- 
nal from diffuse neutral hydrogen towards the end of reion- 
ization. Provided that these peaks are at least marginally 
resolved, the statistical anisotropy of their shapes can be 
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Figure 4. Ratios of the moments, ai/ipi (I = 0, 2, 4), as a function 
of separation. Also shown is the moment aa multiplied by a factor 
of 10. All curves correspond to the linear correlation function of 
ACDM with Q m = 0.3, assuming that the correlation of the mask 
is unity. 

used to constrain the cosmological parameters. This appli- 
cation is along the lines of the classic AP test (Alcock & 
Paczyhski 1979). Since the peaks are compact, foreground 
contamination is minimized in this way of applying the test. 
The velocity flow around the peaks is in the mildly non- 
linear regime and non-linear methods for modeling redshift 
distortions are needed. 

The correlation function of the discrete distribution of 
the peaks could be computed. The anisotropy of this corre- 
lation could be used as a probe of geometric distortions (e.g. 
Ballinger, Peacock & Heavens 1996; da Angela, Outram & 
Shanks 2005) . Here again, foreground contaminations are ir- 
relevant since the 21 cm maps are used merely for locating 
the peaks. One can also consider neutral gas in minihalos in 
which collisions could be efficient at raising the spin temper- 
ature above that of the CMB (Iliev et al 2002) . They amount 
to significant 21 cm emission even before reionization. Here 
also once can apply the AP test to the distribution of these 
minihalos in 21 cm maps. The distribution of these peaks 
would also be affected by redshift distortions. The signature 
of these distortions here is also distinguishable from that of 
geometric distortions. 



6.2 The patchy mask: ionization fronts 

In a UV radiation dominated reionization, the mask made 
by ionized regions is patchy (see figure 3 of Ciardi & Madau 
2002). The boundaries of these patches, i.e., the ionization 
fronts are marked by a sharp transition in 21 cm emission 
and should be possible to identify in future 21 cm maps. 



The shape of the ionization fronts could be used to probe 
geometric distortions. As we have seen in §5.1, the effect 
of peculiar motions on the anisotropy of the boundaries is 
negligible for large enough bubbles. Methods for identify- 
ing these boundaries in noisy maps observed with a given 
instrumental resolution must be developed. This is worth- 
while since foreground contamination are not likely to affect 
the identification of these boundaries and also, by restricting 
the application to large ionized regions, redshift distortions 
can be made negligible. 



7 DISCUSSION 

We have discussed the application of the Alcock-Paczyhski 
test to maps of the redshifted 21-cm emission from the era of 
reionization. Our two main conclusions are a) at z ~ 20 the 
ratio of the frequency to angular distance scales is sensitive 
to the assumed cosmological parameters and the equation of 
state of the dark energy, and b) that geometric distortions 
resulting from an erroneous choice for the cosmological pa- 
rameters are very distinct from redshift distortions caused 
by peculiar motions. Redshift distortions enhance the clus- 
tering along the line of sight while geometric distortions sim- 
ply compress the density contours either along line of sight 
or the perpendicular direction. Judging by the ratio of the 
moments (see figure 4) a determination of the correlation of 
temperature fluctuations to an accuracy of 20% should al- 
low a successful application of the AP test. Prior to the end 
of reionization by UV sources large unionized regions should 
be easily distinguishable from the propagating ionized bub- 
bles in 21 cm emission maps. Temperature correlation func- 
tions over scales up to several Mpcs can be computed re- 
liably from these regions alone. Thus an application of the 
AP test to these correlations completely avoid dealing with 
the mask, M. Nevertheless, for the sake of completeness we 
have presented results including the mask. The prospects for 
measuring the cosmological 21-cm signal are good in view of 
the planed telescopes like the Low Frequency Array t (LO- 
FAR), the Primeval Structure Telescope § (Peterson, Pen & 
We 2005), and the Square Kilometer Array T It remains to 
be seen how well the correlation function will be determined 
from data provided by these telescopes. 

A potential obstacle in extracting any cosmological sig- 
nal from future 21 cm maps is contamination by bright ex- 
tragalactic sources and Galactic synchrotron and free-free 
emision. (e.g. Oh & Mack 2003; Di Matteo et al. 2002; Di 
Matteo, Ciardi & Miniati 2004, Wang et al. 2005). Galac- 
tic contamination contributes about 70% (e.g. Wang et al. 
2005). This component is believed to be a smooth function of 
frequency. If so then it would be possible to remove its con- 
tribution to the correlation function (or line of sight power 
spectrum) of the cosmological signal on small scales. The 
extragalactic point source contamination is more fluctuat- 
ing than the Galactic component. But, as has been shown 

t http://www.lofar.org 

§ http://astrophysics.phys.cmu.edu/ jbp 

^ http://www.skatelescope.org 
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by Wang etal (2005), the total foreground contamination can 
be almost completely cleaned from the line of sight power 
spectrum of the 21 cm emission on scales with wavenumbers 
<C O.lh/Mpc. Therefore, over a wide range of scales, fore- 
ground contaminations are not expected to pose a major 
problem for the application of the AP test and for extract- 
ing other cosmological information from 21 cm maps. 

We have also outlined alternative method for applica- 
tion of the AP test. A promising method relies on identifying 
high density peaks in 21 cm. These peaks could be associ- 
ated with minihalos which could lead to significant emission 
even before the onset of reionization. They could also result 
from neutral hydrogen in dense regions engulfed by the ex- 
panding ionized bubbles during the late stages reionization. 
Another method relies on the geometric distortions of the 
boundaries of patchy ionized regions, i.e. the reionization 
fronts. The advantage of these method is that they are not 
affected by foreground contamination since the 21 cm maps 
are used merely to locate the position of the peaks in the 
first method and the reionization fronts in the second. 

A proper testing of the methodology developed here is 
beyond the scope of the current paper. Ideally one would like 
to use radiative transfer simulations of the reionization era. 
The problem, however, is that most of the simulation boxes 
are still too small for a proper statistical study of reioniza- 
tion. The largest simulation available to us is that of Ciardi 
& Madau 2003, having a cubic box of 20 h~ Mpc on the 
side. We intend to use this simulation to study the distor- 
tions in the correlations and the shapes of the ionization 
fronts. An alternative to full radiative transfer simulations 
is the approximate methods developed by Benson et. al 2001 
to model reionization using large N-body simulations. These 
methods will enable us to study the proposed application of 
the AP test to simulations of box size of ~ 100 h _1 Mpc. 
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